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Abstract 

We demonstrate that in the presence of mixing between different scalar 
leptoquark multiplets it is possible to simultaneously account for the HERA 
high-Q 2 neutral current anomaly, and produce a charged current anomaly of 
comparable magnitude. The reduced branching ratio to electrons and jets of the 
lightest leptoquark state results in a significant weakening of the CDF /DO limits 
on scalar leptoquarks; masses consistent with the HERA neutral current excess 
are comfortably within the allowed range. We show that the possibilities for 
such a successful mixed leptoquark scenario are quite limited, and we investigate 
some aspects of their phenomenology. 
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1 Introduction 



The HI IJ and ZEUS 0] collaborations at HERA recently announced an anomaly at 
high-Q 2 in the e + p — ► eX neutral current (NC) channel. With a combined luminosity 
of 34.3 pb _1 in e + p — > eX mode at \/s = 300 GeV, the two experiments observe 24 
events with Q 2 > 15000 GeV 2 against a Standard Model (SM) expectation of 13.4 ± 
1.0, and 6 events with Q 2 > 25000 GeV 2 against an expectation of only 1.52 ± 0.18. 
The high-Q 2 events seem to be clustered at Bjorken-x values near 0.4 to 0.5, with the 
HI data showing a more pronounced peak. 

HERA is capable of running in two modes: e~p and e + p. In the former mode 
HI and ZEUS have accumulated 1.53 pb~ of data but have observed no statistically 
significant deviations from the SM. Further, the experiments differentiate between 
final state eX and uX, where the neutrino is detected through its missing px- HI has 
also announced its findings in the e + p — ► vX charged current (CC) channel. They 
find 3 events at Q 2 > 20000 GeV 2 with an expectation of 0.74 ± 0.39, but no events 
with Q 2 > 25000 GeV 2 . ZEUS has not announced its CC data. Compared to the NC 
channel, the present CC signal is considerably less statistically significant. However, if 
the current trend persists in the 1997 HERA data, there will be interesting constraints 
on theoretical interpretations of the HERA excess. 

There are three general categories of explanation for the NC excess: (1) a statisti- 
cal fluctuation, (2) unexpected SM physics, for example, a modification of the parton 
distribution functions of the proton at moderate-to-high- x and large Q 2 , or (3) new 
physics. 

For case (1) we have nothing to say. The possibility of modifying the parton 
distribution functions (2) was suggested in ||, and its consequences for the high-Q 2 
HERA NC and CC data were studied in M . This suggestion has the advantage that 
it is relatively conservative, and furthermore automatically avoids the serious flavor 
problems ]5|, |], |7| associated with most new physics explanations. However the most 
attractive hypothesis for modifying the parton distribution functions, the (symmetric) 
intrinsic charm scenario, has already been shown to be inconsistent with the high-Q 2 
HERA data itself |§. 

In category (3) are a variety of forms the new physics could take. These include 
contact interactions |], ^|, ||, s-channel production of a leptoquark or a squark with 
.R-parity violating interactions [§, ||, |^, [11], [12| , and related proposals |13| . 

In this paper we consider the implications of the present and future HERA CC data 
for new physics. The constraints on, and the consequences of, contact interactions 
for the CC data has been considered by Altarelli, et al | JL4|j ; there they found that 
it was very difficult to explain any excess in the CC comparable to that in the NC. 
They also considered CC signals in scenarios with leptoquarks, with emphasis on the 
supersymmetric i?-parity violating case. For the case of leptoquarks they were able 
to find scenarios with a significant CC excess. In this paper we will also consider 



leptoquarks, but unlike Ref. |T4[, we will examine the consequences of leptoquark 
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Table 1: List of scalar leptoquark operators. For each operator, the SU(3) x SU{2) x 
U(l)y quantum numbers of the leptoquark, $, are shown. The fifth column indicates 
whether HERA should find CC events, and the final column lists the mode in which 
HERA should dominantly produce the given leptoquark. Q and L represent SU(2) 
doublet quarks and leptons, while e, u and d are SU (2) singlets. 



mixing on the CC signal and demonstrate that comparable CC and NC excesses can 
similarly be obtained. We also investigate the consequences of leptoquark mixing 
for low-energy precision measurements, and for high-energy direct searches at the 
Tevatron. In this regard, one noteworthy feature of scenarios with a CC signal is that 
the recent CDF/DO leptoquark limits [15, 16, [IT] (for scalar leptoquarks M > 210 GeV 
at 95% C.L.) are considerably weakened. 

If the NC anomaly at HERA is explained by s-channel production of a scalar 
leptoquark then the measured Bjorken-x distribution translates into a mass deter- 
mination: m 2 = xs. For both HI and ZEUS two independent determinations are 
possible, depending on whether x is calculated from the double-angle or electron 
methods. For the seven events selected for special study by HI we find these two 
methods give M 2a = 202 ± 14 GeV and M e = 199 ± 8 GeV respectively. Similarly, 
for the five events selected by ZEUS, M 2a = 231 ± 16 GeV and M e = 219 ± 12 GeV. 
Similarly, the 3 highest-Q 2 events in the HI CC data are also clustered in x yielding 
Mqc — 197 =t 20 GeV, where the error is dominated by the HI energy resolution. 
The fact that the CC and NC mass determinations at HI so closely match supports 
the hypothesis that both anomalies are coming from production of the same particle, 
which then decays either to electrons or neutrinos. 



2 Mixing of scalar leptoquarks 

We will argue in this section that although with a single scalar leptoquark multiplet 
one does not expect a charged-current signal at HERA, in the presence of mixing of 
different leptoquark multiplets a charged-current signal is possible. 

The possible renormalizable couplings of scalar leptoquarks to SM fermions are 
enumerated in Table [I] H], along with their SU (3) x SU (2) x U (l)y quantum numbers. 
We will consider only s-channel production of leptoquarks. (It was shown in |J that 
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w-channel exchange of a light leptoquark could not consistently explain the HERA 
NC data.) If the NC events at HERA are due to an s-channel leptoquark, then the 
scattering of the e + must be off one of the valence quarks in the proton.Q Operators 
satisfying this constraint are listed as "e + " mode in the last column of the table. On 
the other hand, leptoquarks in the remaining "e~" operators are dominantly produced 
in e~ mode, and would have been already detected in HERA's e~p data (even with 
the much smaller accumulated luminosity) given the size of coupling necessary to 
explain the anomaly in the e + p data. The reason for this is that the leptoquarks in 
"e~" operators are produced by scattering off sea quarks in e + p mode, but off valence 
quarks when HERA runs in e~p mode. Finally, in the column "CC" is indicated 
whether or not the operator leads, in the absence of mixing, to CC events at HERA 
in addition to the NC. We see that the three allowed types of scalar leptoquark are 
all SU(2) doublets, and cannot lead to a CC signal at HERA in either e + p or e~p 
mode. 

However, as we now argue, in the presence of mixing, these three operators can 
lead to a CC signal comparable in magnitude to the NC signal. If we wish to simul- 
taneously explain an excess in both the NC and CC channels, then the mixing must 
involve at least one of the following leptoquark states: (i) the charge Q = 2/3 com- 
ponent of $Ld, (ii) the Q = 5/3 component of and (hi) either the Q = —5/3 or 
the Q = —2/3 components of $ e g. These components can mix with other leptoquark 
states or with each other. 

The most important constraint on such mixing arises from the necessity of avoiding 
helicity-unsuppressed rare decays. If a leptoquark couples significantly to both left- 
handed and right-handed leptons and quarks, it generates effective 4- Fermi operators 
of the form u^difiRVi,. These lead, for example, to helicity-unsuppressed n — > ev 
decays, which are severely constrained by experiment [R This eliminates the 
possibility of mixing $ e g with any leptoquark that would lead to a CC signal, since 
by definition such a state has to couple to the left-handed lepton doublet containing 
a v. Similarly, option (ii) is not realized, given that the only other Q = 5/3 state is 
in $*q. Thus we are left with only option (i). 

The Q = 2/3 component of can only mix with either the Q = 2/3 component 
of or with the Q = 2/3 component of the triplet &*lq- Such mixings can arise 
from terms in the Lagrangian, 

C D X^LdH 2 + ~^ Ld § LQ H\ (1) 

where H is the SM Higgs doublet. In a supersymmetric context, we can identify 
H* = Hd and H = H u . The second term arises generically in a softly-broken su- 
persymmetric model, for instance from A-terms. The first term seems to require 
additional fields. 

-'^More specifically, we are excluding scattering off of u or d quarks. We are not considering the 
possibility of scattering off of s or c sea quarks |Io| . ^, where the required Yukawa couplings 
must be of 0(1). 
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In Eq. ([I]) the first term is B and L conserving, while the second term violates L 
(but not B) by AL = 2. Thus the second term can lead to lepton number violating 
processes such as n + — ► e + V e which is not helicity-suppressed. The constraints on the 
AL = process 7r + — > e + v e also apply here, requiring either leptoquark masses over 
several TeV or Yukawa couplings of §lq too small to be of use for producing a CC 
interaction at HERA. Thus we will not consider this possibility any further. 

This leaves us with only the possibility of (<&l<2, §l u ) mixing. In principle both 
the light and heavy mass eigenstate admixtures of the Q = 2/3 leptoquark states 
could be currently produced by HERA. This is because both have couplings to e + d 
after mixing, and it is possible that they have comparable masses. It is also possible 
to produce the <p 5 /^ state at HERA which leads to no CC signal on its own. (We 
will return to the constraints on such a scenario arising from CDF/DO later.) We 
think this is an interesting possibility, but we will focus our discussion on the simpler 
possibility that only the lighter mixed eigenstate is being produced at HERA; all the 

2 /3 

other states are assumed to be heavy. We denote this light mixed state as 4>i = 
cos6<j) 2 /d + sm9(j) 2 /^ , and the heavier mixed state (the orthogonal combination) as 

2 /3 

4>2 . The initial leptoquark Yukawa couplings before mixing are 

£ = \ 1 LM Ld + \ 2 Lm Lu + h.c. (2) 

This leads to couplings of the light (j>\ state of (XicosOed + A 2 sin Bvu)$y Z , and an 
associated CC branching ratio 

Br - (01 ^ ) = A?cos^ + A|sin^ (3) 

For A2 sin# ~ Ai cos# this results in comparable CC and NC signals at HERA. Note 
that to maintain the overall NC event rate, the effective coupling Ai cos 9 (which was 



~ 0.04 in the absence of mixing Q) has to be increased by a factor of 1/J (1 — Br^j). 

It should be noted that CDF has recently announced [0 a lower limit of 210 GeV 
(at 95% C.L.) on scalar leptoquarks that decay into electrons and quarks with branch- 
ing ratio 1, using the NLO QCD calculations of Ref. |17| . (The limits on a leptoquark 
decaying to a quark and a neutrino are much weaker.) This is starting to significantly 
constrain the (unmixed) leptoquark explanation of the HERA NC data. (From the 
calculations of Ref. , it can be seen that this translates into a limit on vector lep- 



toquarks that essentially excludes this option for explaining the HERA NC anomaly.) 

In the mixed scalar leptoquark scenario the situation is significantly changed, 
since the branching ratio to electrons and jets (Br e j = 1 — Br^) is less than one. 
For example, if the branching ratio is 1/2, we estimate the mass limit on to be 
somewhat below 190 GeV, which is perfectly compatible with the HERA data. In 
this case, however, the heavier a ^ so decays into e + + jets with a branching ratio 



2 



This includes the enhancement coming from NLO QCD effects [n2l E(J 
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of 1/2, so it contributes to any possible signal at CDF/DO. Setting the mass of 



<2/3 



at 200 GeV, we estimate that the mass of tfi 2 / 3 must be > 225 GeV (for Br ej - = 1/2) 
using the analysis of Jl7| . 

One consequence of the CDF/DO data is that it strongly disfavors $ 6 q as a possible 
explanation of the HERA NC anomaly As we argued above, it is not possible to 
substantially mix this state with any other which has couplings to neutrinos. Thus 
both components of $ e Q always decay into electrons + jets. Furthermore, since $ e g 
is an unmixed doublet the p-parameter constraint does not allow the mass of the 
heavier component to be arbitrarily large. Typically, one finds that the heavier mass 
is constrained to be below 250 GeV, given that the lighter mass is fixed at ~ 200 GeV. 
Furthermore, there is a perturbativity constraint on the splitting since it is entirely 
due to electroweak symmetry breaking: 5m 2 ~ Xv 2 < (250 GeV) 2 for A < 1. Therefore 
CDF/DO would see two states each decaying to electrons with branching ratio one. 
For instance a heavy state at 250 GeV contributes 0.03 pb to the CDF cross section, 
which implies that the lighter state must be above 220 GeV at 95% C.L. Alternatively 
the quoted CDF limit for a single state implies a limit on a degenerate doublet of 
M$ > 234 GeV at 95% C.L., probably beyond the interesting mass range for HERA. 



2.1 The /^-parameter 

It is natural to consider the implications of the p-parameter constraints on scenarios in 
which scalar leptoquark multiplets mix. We have investigated this in detail following 



Ref. pTJ. In particular Ap = (n + _ — H^/m^r where the indices label SU(2) weak- 



isospin, and the general expression for the quantities Ii a b is given by 

= / |^!Tr([T a , A(k)][T b , A(k)]). (4) 

In this equation the T's are SU(2) matrices in the (reducible) representation of the 
scalar multiplets, A(k) = l/(k% + M 2 ), with M 2 the mass squared matrix, g is the 
weak coupling, and C is a possible color factor (C = 3 for scalar leptoquarks). 

The resulting analytic expressions for Ap in the presence of mixing between two 
SU(2) representations are quite complicated. However, in the case where two weak 
doublet leptoquarks (such as &l<i and &lu) are mixed in such a way that three of the 
four resulting leptoquarks are degenerate, the expressions simplify greatly. In this 
special case 

9 9 2m 2 ml , ( m 2 \ 
m 2 l + m 2 + 2 In — i , 5 

mf — m\ Vml / 

and it is clear then that for maximal mixing (6 = 7r/4), Ap vanishes. For the general 
case we have investigated the contribution of Eq. (f|) to Ap numerically. In many 
cases it is possible to find substantial negative contributions to p. For example, if 
we take the mass eigenvalues to be (202, 273, 240, 240) GeV, with maximal mixing 



3G F cos 2 2fl 

Ap = 



8y/2ir 
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between the first two (the Q = 2/3 states), then we find Ap = —0.001. It is possible 
to find negative values of greater magnitude (especially for large mixing), as well as 
small positive values, if one allows the mass spectrum to vary consistent with the 
CDF/DO limits. Thus the p-parameter provides no significant constraint on the mass 
spectrum in the ($Ld, $l«) mixed case. (Note that if we require substantial mixing 
and one eigenvalue fixed around 200 GeV, then perturbativity does not allow any of 
the other states to become arbitrarily massive.) 



2.2 Atomic parity violation 

The mixing scenarios also have consequences for the atomic parity violation (APV) 
experiments. In the presence of scalar leptoquarks the predicted weak charge of a 
nucleus shifts by 

For the leptoquarks of interest, the values of (5z, #jv) are: (2,1) for and (1,2) for 

Recently the measurement of the weak charge Qw of Cesium (Z = 55, N = 78) 
has significantly improved leading to a new value Q^ pt (Cs) = -72.11 ± 0.93 



This is to be compared to the SM prediction of Q^(Cs) = —73.20 ± 0.09. There 
are also reasonably good measurements for Thallium {Z — 81, N — 205), leading to 
Q^ pt (Th) = -115.0 ±4.5 compared to the prediction Q^(Th) = -116.8 ±0.19 [[23j. 
For the central values of the masses and couplings of the and <&Ld states found 
in a fit to the NC HERA data (X Lu ~ 0.020, X Ld ~ 0.04, with masses ~ 200 GeV in 
both cases), Qw(Cs) shifts by the small amounts AQ^(Cs) = —0.09 and AQ^(Cs) = 
—0.40, respectively, further away from experiment. 

In the presence of mixing, the expression for AQw changes to 



(g w /M w ] 
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'AiCos#\ 2 /Aisin6>\ 2 
mi j \ m 2 J 



+ ■■■ (7) 



The ellipsis denotes the contributions from the unmixed components of the heavier 
multiplet, which depend upon its unknown Yukawa coupling A2 to SM fermions. Since 
in the ($Ld, &Lu) mixing case the combination Ai cos 9 is now fixed to be ~ 0.04 by 
the magnitude of the NC excess, it is clear that the (Ai sin#/m2) 2 term increases the 
shift in AQw compared to the unmixed case. For instance, in the case of maximal 
mixing one finds roughly a doubling of AQw, relative to the unmixed case. (This shift 
is only increased by the terms proportional to \\ omitted in Eqn. (|7]).) For Thallium, 
all the shifts are about twice as large, but because of the much larger experimental 
error this corresponds to a statistically less significant change. 

Thus, overall, the discrepancy between the APV measurements and the SM tends 
to increase in the case of significant ($Ld, $l«) mixing, by approximately 0.75cr for Ce- 
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sium (the exact value depending on the mass the mixing angle, and the unknown 
Yukawa coupling A 2 ). 



2.3 Heavy leptoquark decays 

We have seen from the discussions above that the leptoquark interpretation of the 
HERA NC data admits a comparably large cross section in the charged current mode 
as well, provided that the leptoquark produced is an admixture of different lepto- 
quark states. The masses of these additional leptoquarks are constrained by the 
p-parameter, and perturbativity, to be roughly below 300 GeV. It is then of interest 
to ask how these extra leptoquarks could be detected once they are produced at col- 
liders. Here we wish to point out some interesting decay channels of these heavier 
leptoquarks which will be of relevance for their search at the Tevatron. 

For (&Ld, &Lu) leptoquark mixing, since both are doublets, there are four physical 
states. The 0^ 3 mixes with the (J) 2 / 3 forming two mass eigenstates (0 2 ^ 3 and 0^ 3 ) 
of masses mi and 1712, with mi f=s 200 GeV. The two other states (f)^ 3 and 4> 5 / 3 
are unmixed, and we denote their masses by m 3 and m 4 . The two body decay 

l jo 

^Ld d + V is always open, but the width for this decay is suppressed by a small 
Yukawa coupling squared A 2 ~ (0.04) 2 . If m 3 is larger than m x + m w , the decay 
<fi L l J 3 — > 2//3 + W will be open and will be dominant. One signature will be 3j + e; 
a more spectacular signature into Ivej occurs about 1/3 as often. If this channel is 
closed, there is still the three-body decay through a virtual W which can compete 
with the d + V decay rate. When the W is off-shell, the rate for the three-body decay 
is given by 

iW - 0/ //) = ^ 3 m 3 I dt — (8) 

where r = m 2 /m|, s = m^/m^, A s = Af = g/y/2, and summation of different 
helicities is assumed. The competing decay has a rate 

r(0 L ] /3 ^d + V) = ^-m 3 . (9) 

Taking mi = 200 GeV, the total width for the three body decay is found to be 
(0.73 MeV, 2.0 MeV, 5.3 MeV) for m 3 = (250, 260, 270) GeV. These numbers are 
to be compared with the two-body decay width T « 7 MeV for A = 0.04. Identical 

5 /3 

results hold for the decay of the Ltt leptoquark. 

What about the heavier (j)^ 3 leptoquark? Since the mass eigenstates are combina- 
tions of Q =2/3 states carrying different T 3 , there is a residual off-diagonal 2//3 c/) 2 / 3 'Z 
coupling with a coefficient [g sin 29] /(2 cos 0w), where 9 is the leptoquark mixing an- 
gle. This leads to the three body decay (f)^ 3 — > 2 ^ 3 + Z* — > (f)^ 3 + f + f. The rate is 
given by Eq. (1), now with A s = gsin29/(2cos9 w ), Af = g(T 3 — Qsm 2 9w)/ costf^, 
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r = m\/m\ and s = m 2 z /m\. Again, this mode can compete with the 2-body decay 
mode. For instance, if m-i = (260, 270, 280) GeV with maximal mixing then the total 
width for the 3-body decay is (0.7, 1.8, 4.1) MeV. If enough of the heavier leptoquarks 
can be produced, the gold-plated signature will be £ + £~e + e~jj. Here the £ + £~ come 
from the decay of 02 to Z*<pi with the 0! then further decaying to ej; the second 02 
on the other hand decays directly to ej. 

3 Conclusions 

In the presence of mixing between different scalar leptoquark SU(2) multiplets we have 
demonstrated that it is possible to simultaneously account for the HERA high-Q 2 NC 
anomaly, and produce a CC anomaly of comparable magnitude. The branching ratio 
to electrons and jet(s) of the lightest leptoquark state that by assumption is being 
produced at HERA is consequently reduced. This results in a significant weakening of 
the recently announced CDF/DO limits on scalar leptoquarks, and masses consistent 
with an explanation of HERA are comfortably within the allowed range. 

However we showed that there is only one possibility for such a successful mixed 
leptoquark scenario: the §Ld doublet leptoquark multiplet must mix with the $l u 
doublet. We investigated various aspects of the associated phenomenology. In par- 
ticular, this scenario often allows relatively large splittings among the leptoquarks 
while satisfying the p-parameter constraints. Indeed it is possible to find substantial 
negative contributions to p in the case of large mixing. We also point out that there 
are novel possibilities for the decays of heavier members of the leptoquark multiplets, 
which might be relevant for searches at the Tevatron. 

Finally, it should be emphasized that the mixed scenario does not address what 
is in our opinion the most serious problem facing leptoquark (or .R-parity violating 
squark) explanations of the HERA anomaly, namely, the very serious flavor problem 
that all such models engender || || [/J. 
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